Multiple firing coherence resonances in excitatory and inhibitory coupled neurons 
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The impact of inhibitory and excitatory synapses in delay-coupled Hodgkin-Huxley neurons that are driven 
by noise is studied. If both synaptic types are used for coupling, appropriately tuned delays in the inhibition 
feedback induce multiple firing coherence resonances at sufficiently strong coupling strengths, thus giving rise 
to tongues of coherency in the corresponding delay-strength parameter plane. If only inhibitory synapses are 
used, however, appropriately tuned delays also give rise to multiresonant responses, yet the successive delays 
warranting an optimal coherence of excitations obey different relations with regards to the inherent time scales 
of neuronal dynamics. This leads to denser coherence resonance patterns in the delay-strength parameter plane. 
The robustness of these findings to the introduction of delay in the excitatory feedback, to noise, and to the 
number of coupled neurons is determined. Mechanisms underlying our observations are revealed, and it is 
suggested that the regularity of spiking across neuronal networks can be optimized in an unexpectedly rich 
variety of ways, depending on the type of coupling and the duration of delays. 
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I. INTRODUCTION 

Neurophysiologies studies have revealed the existence of 
accurately timed patterns of spikes by a variety of cogni- 
tive and motoric tasks The timing of these spikes, or 
neuronal firings, is accurate to within the millisecond range, 
which poses great challenges with regards to the identification 
of mechanisms that would be able to ensure such precision. 
Following their initial observation in the cortex of monkeys 
UIH], me precisely timed spikes have been reported and in- 
vestigated for motor functions [3], the neuronal response of 
visual systems [4], and the complex spatial fingertip events 
to name but a few examples. Not surprisingly, synchro- 
nized, precisely timed firings can be observed at virtually all 
neuronal processing levels, including the retina [9], the lateral 
geniculate nucleus IToll . and the cortex iflll, [l2tl . 

Since it is well known that noise can play a constructive 
role in different types of nonlinear dynamical systems, which 
arguably describe also neuronal dynamics 11311 . this opens the 
possibility of exploiting such mechanisms for explaining, or 
at least supporting, the aforementioned precision of neuronal 
firings. Stochastic resonance Il4l4l6ll and coherence reso- 
nance lfl7l - tl9ll are amongst the most prominent examples by 
means of which noise of appropriate intensity is able either 
to enhance the detection of weak deterministic signals [20] or 
evoke coherent response in nonlinear dynamical systems in 
the absence of any deterministic inputs. The potential benefits 
of noise range from ice ages to crayfish and SQUIDs [21], to 
neural systems, as most recently reviewed in [22]. 

Following initial advances on individual dynamical sys- 
tems, the focus begun shifting to spatially extended systems 
l23tl . especially also to such with complex networks describ- 
ing connections between the individual units J24l l25ll . For 
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example, coherence resonance on a small world network was 
investigated in lErjll . while array-enhanced resonances were 
reported in lf27ll . Moreover, spatial coherence resonance was 
observed first near pattern-forming instabilities II28I1 . and lat- 
ter also in excitable media [29]. Excitable systems in general 
proved to be very susceptible to a multitude of noise-induced 
phenomena, as reviewed comprehensively in (30J. Adding 
spatial degrees of freedom, along with the possibilities for in- 
troducing other sources of heterogeneity, lead to the discovery 
of very interesting and quite exotic phenomena, such as the 
ghost resonance 113 ill , and double as well as multiple stochas- 
tic lH^-fHtl and coherence ll36l - l39ll resonances. 

For neural systems, a wealth of interesting and new phe- 
nomena was made observable by integrating realistic features 
of neuronal dynamics into the studied models. Information 
transmission delays or synaptic delays, for example, are inher- 
ent to the nervous system because of the finite speed at which 
action potentials propagate across neuron axons, and due to 
time lapses occurring at both dendritic and synaptic process- 
ing 14011 . Following seminal works examining the impact of 
delays on excitable and other dynamical systems [41-43], the 
stability and attainability of synchronous oscillations ||44 -|46| 
and the role of delays in shaping spatiotemporal dynamics of 
neuronal activity [47] were investigated. Moreover, the role of 
delays in coupled Hodgkin-Huxley neurons was also investi- 
gated for the phenomenon of coherence resonance, and it was 
reported that properly tuned delays can lead to the occurrence 
of multiple resonances j48il49ll . 

In this letter, we extend the scope of coherence resonance in 
models of neuronal dynamics by considering besides synap- 
tic delays also different types of synaptic coupling. While 
the role of chemical synapses in coupled neurons with noise 
has been investigated in [50], and although the general dy- 
namics of sparsely connected networks of excitatory and in- 
hibitory spiking neurons is known |Sl|] , our approach, joining 
these distinctive features of neuronal dynamics (synaptic de- 
lays, different types of synaptic coupling, and noise), allows 
for the identification of new ways by means of which the co- 
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herence, and thus the accuracy of neuronal firings, can be im- 
proved. Most interestingly, we report the occurrence of multi- 
ple coherence resonance patterns in the corresponding delay- 
strength parameter plane when either inhibitory and excitatory 
or only inhibitory synapses are used for coupling. The details 
of these multiple firing coherence resonances, and in particu- 
lar the conditions at which they occur, however, depend sig- 
nificantly on the type of coupling. Reported results suggest 
that characteristic time scales related to the information trans- 
mission and inhibition in neuronal networks may interplay in 
intricate ways, and by doing so give rise to new mechanisms 
for optimizing spiking regularity. 

The remainder of this letter is organized as follows. In the 
next section we describe the model, then we present the main 
results separately for the two coupling scenarios, while lastly 
we summarize our findings and discuss their potential impli- 
cations. 



II. MODEL DEFINITION 

For simplicity, we consider two Hodgkin-Huxley neurons 
lfl3ll that are coupled by inhibitory and/or excitatory synapses. 
Equations describing the dynamics are: 

C ^TT = -9Nam 3 h(V-V Na )-g L (V-V L ) 
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where Vi is the transmembrane potential of the z-th neuron. 
Moreover, rrii, hi and rij are the gating variables, where the 
voltage-dependent opening and closing rates are: 
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The membrane capacity is C = 1 (^F/cm 2 ), and g^ a = 120 
/iF/cm 2 , gx = 36 /iF/cm 2 and g^ = 0.3 /iF/cm 2 are the max- 
imal sodium, potassium and leakage conductances, respec- 
tively. The corresponding reversal potentials are Vjsr a = 50 



mV, Vk = — 77 mV and Vl — —54.4 mV. Using these pa- 
rameter values, a single Hodgkin-Huxley neuron has a sub- 
critical Hopf bifurcation at the external current I = Ii = 
9.8/iA/cm 2 . Between J = J 2 = 6.2/iA/cm 2 and I\ sta- 
ble limit cycles coexist with stable steady states, whereas for 
I < I2 (I > Ii) excitable steady states (limit cycles) are the 
only stable solutions. If / > 155/iA/cm 2 , on the other hand, 
the oscillations vanish by means of a supercritical Hopf bifur- 
cation. A more detailed bifurcation analysis of the Hodgkin- 
Huxley model was performed in lf52l p3tl . Here we are in- 
terested in the region I < I 2 , where neurons are unable to 
fire spontaneously, i.e, remain forever quiescent in the ab- 
sence of external stimuli. We thus set I = 6.1/iA/cm 2 , so 
that both neurons are in an excitable steady state. Gaussian 
noise £i(t), having mean < &(t) > = and autocorrelation 
< £i(t)£j (t') >— SijS(t —If), thus acts as the source of large- 
amplitude excitations, where a determines the noise intensity. 

We consider two different coupling schemes. First, the two 
neurons are coupled in a hybrid way using inhibitory and ex- 
citatory synapses. The coupling terms in this case are: 
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where the inhibitory feedback is delayed by t. Second, only 
inhibitory synapses are used for coupling, in which case the 
coupling becomes: 



p,3 = 

syn 



ginh 



(Vi - V mh ) 



(l+exp{-A[V^-T)-9 s ]}) : 
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where the inhibitory feedback is again delayed by r, only 
that here this applies to both directions. In the above cou- 
pling terms gi n h( eX c) determines the strength of the synaptic 
conductance, i.e., the coupling strength, while Vi n h = —80 
mV and V exc = 20 mV are the reversal potentials for the in- 
hibitory and the excitatory synapse, respectively. Moreover, 
S = is the threshold, above which the postsynaptic neuron 
is affected by the presynaptic one, and A = 10 is a constant 
rate for the onset of excitation or inhibition. In what follows, 
we will investigate the impact of the delay r and the coupling 
strength gmhfexc) on me occurrence of firing coherence res- 
onance, and we will do so separately for the two described 
coupling schemes. 



III. RESULTS 

We start by presenting the results as obtained with hybrid 
coupling, i.e., when excitatory and inhibitory synapses are 
used for connecting the two neurons. Figure [1] features char- 
acteristic time courses of the transmembrane potential V of 
the excitatory neuron, from where it can be observed at a 
glance that the coherence of excitations depends critically on 
the delay of the inhibitory feedback r. Importantly though, 
the relation between the coherency and the value of r is not 
monotonous, but rather it is intermittent. That is to say, as r 



3 




600 800 1000 1200 1400 

f 

FIG. 1: Appropriately adjusted delays r in the one-directional in- 
hibition feedback enhance the regularity of spiking by hybrid cou- 
pling of the two neurons. Depicted are characteristic time courses of 
the transmembrane potential V of the excitatory neuron for different 
values of r: (a) 0, (b) 8.0, (c) 20, (d) 24, (e) 35 and (f) 40. It can 
be observed that the regularity of spiking in panels (b), (d) and (f) 
(traces depicted green) is higher than in panels (a), (c) and (e) (traces 
depicted red). Other parameter values are: g exc = 0.11, gi n h = 1.0 
and a = 1.5. 



increases the regularity is lost and regained intermittently as 
different values of r come to determine the delay of inhibi- 
tion. Time courses depicted green (panels b, d and f) exhibit 
more coherent spiking than time courses depicted red (pan- 
els a, c and e). This is characteristic for multiresonant phe- 
nomena, and in fact these observations can be made quanti- 
tatively more precise by introducing a coherence measure C 
as follows. Let the sequence to < t\ < < ■ ■ ■ < t n 
denote the firing times of the considered neuron. From the se- 
quence of {tk}, the interspike intervals (ISI) are determined 
as Tfc = tk — tk-i(k = 1, 2, • • • , n). To characterize the co- 
herence of the firings, the measure C is defined as 

C= ^< T Z>-< T *>1 (14) 
<T k > 

where (•) is the time average. In particular, C is the ratio of 
the standard deviation and the average of the interspike in- 
tervals, and it is indeed an excellent quantity for effectively 
determining the occurrence of coherence resonance from neu- 
ronal firing. From Eq. (14) it follows that the more coherent 
the firing, the smaller the value of C. We would also like to 
note that C is the reciprocal of the coefficient of variation in 



FIG. 2: Delay-induced multiresonances in case of hybrid coupling 
of the two neurons. Panel (a) shows the coherence measure C in 
dependence on r for different values of gi n h- It can be observed 
that the stronger the coupling the better expressed the recurrently ap- 
pearing minima of C. Panel (b) features the contours of C (white 
depicts minimal and black maximal values) on the corresponding 
delay-strength r — g in h parameter plane, where multiple tongues of 
coherency (white) emerge due to an interplay between the synaptic 
delay r and the characteristic time scale of the two Hodgkin-Huxley 
neurons (as determined by the characteristic excitatory time T e and 
the complex conjugate part of the eigenvalues of the excitatory steady 
state). Other parameter values are: a = 1.5. 



a point process, which is widely used in the field of neuro- 
science 15411 . 

Using the introduced coherence measure C, we demon- 
strate in Fig. |2] the occurrence of multiresonant behavior in 
dependence on r. Results presented in panel (a) indicate that 
C has several minima in the considered interval of r, and that 
these are better pronounced, i.e., less susceptible to statisti- 
cal deviations, for larger coupling strengths ginh- In general, 
however, the dependence of C on gi n h is fairly insignificant, 
pointing towards the fact that in case of hybrid coupling the 
strength of the synaptic conductance of one type (e.g., the in- 
hibitory type) has little impact if the other (e.g., the excitatory 
type) remains unchanged. The contours in panel (b) confirm 
this, as the tongues of coherency (white regions) simply shrink 
in width as <?j n /j decreases, but otherwise do not alter the de- 
pendence of C on the inhibition delay t. In many ways, these 
results are reminiscent of delay-induced multiple stochastic 
resonances that were previously reported for scale-free neu- 
ronal networks [33], and are indicative for an interplay be- 
tween the time scales inherent to the system dynamics and the 
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FIG. 3: Appropriately adjusted delays r in the bidirectional inhibi- 
tion feedback enhance the regularity of spiking by inhibitory cou- 
pling of the two neurons. Depicted are characteristic time courses of 
the transmembrane potential V of one neuron for different values of 
r: (a) 0, (b) 2.0, (c) 5.0, (d) 11, (e) 15 and (f) 19. As in Fig.Q it 
can be observed that the regularity of spiking in panels (b), (d) and 
(f) (traces depicted green) is higher than in panels (a), (c) and (e) 
(traces depicted red). Other parameter values are: g in h — 0.75 and 
a = 1.5. 



time scales introduced by means of the delay. 

Turning to the second coupling scheme relying only on in- 
hibitory synapses, however, we find somewhat unexpected re- 
sults. While the time courses of the transmembrane potential 
V presented in Fig. [3] do not suggest quantitatively different 
behavior in that certain values of r warrant higher coherency 
of spiking than other values (which is also what we can ob- 
serve in Fig. [TJ, a more accurate quantitative analysis pre- 
sented in Fig. |4] indicates otherwise. In particular, in panels 
(a) and (b) we find that the minima of C are much more fre- 
quent in the considered span of t values as this was the case 
for hybrid coupling. While for the later a total of three minima 
can be observed within < r < 50 (see Fig. [2j, for purely in- 
hibitory coupling twice as many minima are inferable within 
the same span of r values. 

The origins of these multiresonant phenomena can be 
linked to different inherent properties of neuronal dynamics. 
First, it is useful to define the so called average excitatory 
time T e , which is the average time between two consecutive 
spikes. For an isolated Hodgkin-Huxley neuron driven by 
noise this time decreases and saturates towards T e w 16 for 
a > 4.0 (note that this corresponds to a strong noise limit, 
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FIG. 4: Delay-induced multiresonances in case of inhibitory cou- 
pling of the two neurons. Panel (a) shows the coherence measure C 
in dependence on r for different values of gi n h- As in Fig.[2] it holds 
that the stronger the coupling the better expressed the recurrently ap- 
pearing minima of C. However, in the considered span of r values, 
twice as many minima as by hybrid coupling can be observed. Panel 
(b) features the contours of C (white depicts minimal and black max- 
imal values) on the corresponding delay-strength r — g in h parameter 
plane, where the much denser tongues of coherency are clearly in- 
ferable. This indicates that the interplay between the synaptic delay 
r and the characteristic time scale of the two Hodgkin-Huxley neu- 
rons is more efficient by purely inhibitory coupling. Other parameter 
values are: a — 1.5. 



above which the system may already exhibit numerical insta- 
bility). Increasing the noise intensity further and lowering the 
time step for numerical integration, it is in principle possi- 
ble to arrive at even lower average excitatory times T e k, 12, 
which agrees with the theoretical prediction stemming from 
the imaginary parts of the complex conjugate eigenvalues 
ImXij = ±icj = ±i0.54, where T e = 2tt/lj = 11.63. Since 
in our simulations, however, we use a comparatively low noise 
intensity a = 1.5, the average excitatory time T e ss 16 of an 
isolated Hodgkin-Huxley neuron is the more accurate approx- 
imation for the inherent time scale of the considered neuronal 
dynamics. For hybrid coupling, we thus find the first mini- 
mum of C at T e /2, and subsequent minima at odd multiples 
of the half of the average excitatory time [see Fig.[2fa)], which 
agrees with the doubly effect of the two considered synaptic 
types l55ll . The average excitatory time is reflected also in the 
time courses presented in Figs. [TJb,d,f) (note that in these the 
firing is accurate and ordered due to the constructive impact of 
t), where the average spiking period is approximately equal 
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FIG. 5: Delay-induced multiresonances in the presence of additional 
delay in the excitatory feedback and noise. Panel (a) features a com- 
parison of the coherence measure C as obtained with and without 
excitatory synaptic delay in dependence on r for hybrid coupling. 
It can be observed that the introduction of delays in the excitatory 
feedback can substantially reduce delays warranting the most coher- 
ent response. Other parameter values are: g exc = 1.0, a — 1.5. 
Panel (b) depicts C in dependence on r for different values of the 
noise intensity a in two purely inhibitory coupled neurons. It can 
be observed that as the intensity of noise increases the maximally at- 
tainable values of C decrease (yet the effect saturates for higher a). 
Optimal delays, however, remain unaffected by noise, which indi- 
cates robustness of the observed delay-induced multiresonances. 



toT e . 

Conversely, for inhibitory coupling, the matching of the 
time scales leading to the multiresonant dependence of C on 
t is different. Although the average excitatory time T e w 16 
is likewise [as in Figs. 01b,d,f)] reflected in the correspond- 
ing time courses presented in Figs. [3|b,d,f), which have the 
same average inter-spike interval, twice as many minima im- 
ply that the resonant matching occurs not just for odd multi- 
ples of T e /2, but in fact for odd and even multiples. However, 
all the minima of C are preceded by a small delay of 2s (where 
the first minimum occurs) that is necessary for the first reso- 
nant response. Since the purely inhibitory type of synaptic 
coupling lacks the excitatory input that is present by hybrid 
coupling, in the former case the matching of the time scales is 
twice as efficient. 

Finally it is of interest to examine the robustness of our 
findings in the presence of delayed excitatory feedback, dif- 
ferent levels of noise, and for different sizes of the network. In 
Fig- 13 a), we present the results with and without delayed ex- 
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FIG. 6: Delay-induced multiresonances in a ring network consisting 
of 100 neurons. Panel (a) features results as obtained with delay in 
the excitatory feedback (g e xc) and hybrid coupling. Panel (b), on the 
other hand, depicts C in dependence on r as obtained with delay in 
the inhibitory feedback (ginh) and purely inhibitory coupling. Based 
on the presented results, it can be concluded that multiresonances in 
a ring network can be observed irrespective of the coupling and delay 
type, if only the delays are appropriately adjusted. However, delays 
warranting optimal coherence in the network with purely inhibitory 
coupling (b) are smaller that those in the network with hybrid cou- 
pling (a). Other parameter values are: g in h = 1-5 [applicable for 
panel (a) only] and a — 1.5. 



citatory feedback in a hybridly coupled two-neuron system. It 
can be observed that, while multiresonances can be observed 
in both cases, the introduction of delays also in the excita- 
tory feedback (in addition to delays in the inhibitory feedback) 
may substantially reduce the delays that warrant an optimal 
response of the system (maximal values of C). Thus, delayed 
excitatory feedback does affect the results quantitatively, yet it 
does not affect the qualitative picture. Figure [2b) shows that 
different noise intensities a have a similar impact. In particu- 
lar, while higher values of a may reduce maximally attainable 
values of C, the multiple maxima are always clearly inferable 
and their positions do not shift. Hence, noise is also unable 
to significantly affect the results. Lastly, we present in Fig. [6] 
results obtained on a larger ring network for the two different 
coupling types. Regardless of whether the coupling is hybrid 
with delays introduced to both types of synapses [panel (a)] 
or purely inhibitory [panel (b)], the multiple coherence reso- 
nances are clearly inferable. Importantly, also on larger net- 
works the purely inhibitory mode of interneuronal communi- 
cation appears to be more efficient (there are more maxima 
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of C in a given span of r) than the hybrid mode, which fully 
agrees with our conclusions obtained by means of the analysis 
of the two-neuron system, and thus solidifies the high robust- 
ness of our main conclusions, which we will summarize in 
what follows. 



IV. SUMMARY 

Summarizing, we have demonstrated the occurrence of 
multiresonant elevation of firing precision, as quantified 
by means of a coherence measure, in synaptically coupled 
Hodgkin-Huxley neurons. We have separately considered hy- 
brid and purely inhibitory coupling, and we have discovered 
that the resonant matching of the different time scales that 
are inherent to the Hodgkin-Huxley model (and the informa- 
tion transmission delay) is twice as efficient in the latter case. 
Our results thus reveal unexpected possibilities for the reso- 
nant enhancement of firing precision by means of matching 
of different time scales of neuronal dynamics. Moreover, we 
have examined the robustness of our findings to the introduc- 
tion of delay in the excitatory feedback, to noise, and to the 



number of coupled neurons. We have found that delayed exci- 
tatory feedback may substantially reduce the length of delays 
that ensure an optimal response of the system, yet that it does 
not qualitatively affect the results. Neither do noise and the 
size of the network, which led us to the conclusion that the 
reported results are highly robust, and that they are thus ex- 
pected to remain valid also in other related neuronal systems. 
We hope that our study will prove useful for facilitating the de- 
velopment of concepts such as function-follow-form ll56ll57tl 
and the application of methods of statistical phy sics for better 
understanding conditions such as epilepsy [58-60] and other 
neurodegenerative diseases, as well as for better understand- 
ing the mechanisms behind high-precision firing patterns in 
more realistic neuronal networks. 



Acknowledgments 

This research was supported by the National Science Foun- 
dation of China (grants 1 1 172017 and 10832006) and by the 
Slovenian Research Agency (grant Jl-4055). 



[1] R. Lestienne, B. Strehler, Time structure and stimulus depen- 
dence of precisely replicating patterns present in monkey corti- 
cal neuronal spike trains, Brain Res. 437 (1987) 214. 

[2] M. Abeles, H. Bergman, F. Margalit, E. Vaadia, Spatiotempo- 
ral firing patterns in the frontal cortex of behaving monkeys, J. 
Neurophysiol. 70 (1993) 1629. 

[3] A. Riehle, S. Grim, M. Diesmann, A. Aertsen, Spike synchro- 
nization and rate modulation differentially involved in motor 
function, Science 278 (1997) 1950. 

[4] M. Oram, M. Wiener, R. Lestienne, B. Richmond, Stochastic 
nature of precisely timed spike patterns in visual system neu- 
ronal responses, J. Neurophysiol. 81 (1999) 3021. 

[5] R. S. Johansson, I. Birznieks, First spikes in ensembles of hu- 
man tactile afferents code complex spatial fingertip events, Nat. 
Neurosci. 7 (2004) 170. 

[6] G. Pipa, A. Riehle, S. Grim, Validation of task-related excess 
of spike coincidences based on neuroxidence, Neurocomputing 
70 (2007) 2064. 

[7] H. J. Cao, Miguel A. F. Sanjuan, A mechanism for elliptic-like 
bursting and synchronization of bursts in a map-based neuron 
network, Cognitive Processing 10 (2009) 23. 

[8] B. Ibarz, Jose M. Casado, A. F. Sanjuan, Map-based models in 
neuronal dynamics, Physics Reports, 501(1-2) (2011) 1. 

[9] S. Neuenschwander, W. Singer, Long-range synchronization of 
oscillatory light responses in the cat retina and lateral geniculate 
nucleus, Nature 379 (1996) 728. 
[10] M. Castelo-Branco, S. Neuenschwander, W. Singer, Synchro- 
nization of visual responses between the cortex, lateral genicu- 
late nucleus, and retina in the anesthetized cat, J. Neurosci. 18 
(1998) 6395. 

[11] C. M. Gray, P. Konig, A. K. Engel, W. Singer, Oscillatory 
responses in cat visual cortex exhibit inter-columnar synchro- 
nization which reflects global stimulus properties, Nature 338 
(1989) 334. 

[12] P. Fries, P. R. Roelfsema, A. K. Engel, P. Konig, W. Singer, 



Synchronization of oscillatory responses in visual cortex cor- 
relates with perception in interocular rivalry, Proc. Natl. Acad. 
Sci. USA 94 (1997) 12699. 

[13] A. L. Hodgkin, A. F. Huxley, A quantitative description of 
membrane current and its application to conduction and exci- 
tation in nerve, J. Physiol. 117 (1952) 500. 

[14] R. Benzi, A. Sutera, A. Vulpiani, The mechanism of stochastic 
resonance, J. Phys. A 14 (1981) L453. 

[15] C. Nicolis, G. Nicolis, Stochastic aspects of climatic transitions 
- additive fluctuations, Tellus 33 (1981) 225. 

[16] L. Gammaitoni, P. Hanggi, P. Jung, F. Marchesoni, Stochastic 
resonance, Rev. Mod. Phys. 70 (1998) 223. 

[17] H. Gang, T. Ditzinger, C. Z. Ning, H. Haken, Stochastic res- 
onance without external periodic force, Phys. Rev. Lett. 71 
(1993) 807. 

[18] A. Longtin, Autonomous stochastic resonance in bursting neu- 
rons, Phys. Rev. E 55 (1997) 868. 

[19] A. S. Pikovsky, J. Kurths, Coherence resonance in a noise- 
driven excitable system, Phys. Rev. Lett. 78 (1997) 775. 

[20] P. Hanggi, Stochastic resonance in biology, how noise can en- 
hance detection of weak signals and help improve biological 
information processing, ChemPhysChem 3 (2002) 285. 

[21] K. Wiesenfeld, F. Moss, Stochastic resonance and the bene- 
fits of noise: from ice ages to crayfish and squids, Nature 373 
(1995) 33. 

[22] M. D. McDonnell, L. M. Ward, The benefits of noise in neural 
systems: bridging theory and experiment, Nat. Rev. Neurosci. 
12(2011)415. 

[23] F. Sagues, J. M. Sancho, J. Garcfa-Ojalvo, Spatiotemporal order 
out of noise, Rev. Mod. Phys. 79 (2007) 829. 

[24] R. Albert, A.-L. Barabasi, Statistical mechanics of complex net- 
works, Rev. Mod. Phys. 74 (2002) 47. 

[25] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, D. Hwang, 
Complex networks: Structure and dynamics, Phys. Rep. 424 
(2006) 175. 



7 



[26] O. Kwon, H.-T. Moon, Coherence resonance in small-world 
networks of excitable cells, Phys Lett. A 298 (2002) 319. 

[27] C. Zhou, J. Kurths, B. Hu, Frequency and phase locking 
of noise-sustained oscillations in coupled excitable systems: 
Array-enhanced resonances, Phys. Rev. E 67 (2003) 030101. 

[28] O. Carrillo, M. A. Santos, J. Garcfa-Ojalvo, J. M. Sancho, Spa- 
tial coherence resonance near pattern-forming instabilities, Eu- 
rophys. Lett. 65 (2007) 452. 

[29] M. Perc, Spatial coherence resonance in excitable media, Phys. 
Rev. E 72 (2005) 016207. 

[30] B. Lindner, J. Garcfa-Ojalvo, A. Neiman, L. Schimansky-Geier, 
Effects of noise in excitable systems, Phys. Rep. 392 (2004) 
321. 

[31] P. Balenzuela, J. Garcfa-Ojalvo, E. Manjarrez, L. Martinez, 
C. R. Mirasso, Ghost resonance in a pool of heterogeneous neu- 
rons, BioSystems 89 (2007) 166. 

[32] A. Zaikin, J. Garcfa-Ojalvo, R. Bascones, E. Ullner, J. Kurths, 
Doubly stochastic coherence via noise-induced symmetry in 
bistable neural models, Phys. Rev. Lett. 90 (2003) 030601. 

[33] Q. Wang, M. Perc, Z. Duan, G. Chen, Delay-induced multiple 
stochastic resonances on scale-free neuronal networks, Chaos 
19(2009) 023112. 

[34] C. Gan, M. Perc, Q. Wang, Delay-aided stochastic multireso- 
nances on scale-free fitzhugh-nagumo neuronal networks, Chin. 
Phys. B 19 (2010) 040508. 

[35] Zeng, C. PL, Gong, A. L., Zeng, C. P., Nie, L. R., Stochas- 
tic multi-resonance in an overdamped bistable system with two 
types of modulation signal, Eur. Phys. J. D 62 (201 1) 219. 

[36] Y. Horikawa, Coherence resonance with multiple peaks in 
a coupled fitzhugh-nagumo model, Phys. Rev. E 64 (2001) 
031905. 

[37] T. Kreuz, S. Luccioli, A. Torcini, Double coherence resonance 
in neuron models driven by discrete correlated noise, Phys. Rev. 
Lett. 97 (2006) 238101. 

[38] W. Bao-Hua, L. Qi-Shao, L. Shu- Juan, L. Xiu-Feng, Spatiotem- 
poral multiple coherence resonances and calcium waves in a 
coupled hepatocyte system, Chinese Physics B 18 (2009) 872. 

[39] X. Lin, Y. Gong, L. Wang, Multiple coherence resonance in- 
duced by time-periodic coupling in stochastic Hodgkin-Huxley 
neuronal networks, Chaos 21 (2011) 043109. 

[40] E. R. Kandel, J. H. Schwartz, T. M. Jessell, Principles of Neural 
Science, Elsevier, Amsterdam, 1991. 

[41] F. M. Atay, Distributed delays facilitate amplitude death of cou- 
pled oscillators, Phys. Rev. Lett. 91 (2003) 094101. 

[42] G. C. Sethia, J. Kurths, A. Sen, Coherence resonance in an ex- 
citable system with time delay, Phys. Lett. A 364 (2007) 227. 

[43] G. C. Sethia, A. Sen, F. M. Atay, Clustered chimera states in 
delay-coupled oscillator systems, Phys. Rev. Lett. 100 (2008) 
144102. 

[44] E. Rossoni, Y. Chen, M. Ding, J. Feng, Stability of synchronous 
oscillations in a system of hodgkin-huxley neurons with de- 



layed diffusive and pulsed coupling, Phys. Rev. E 71 (2005) 
061904. 

[45] Q. Wang, Z. Duan, M. Perc, G. Chen, Synchronization tran- 
sitions on small-world neuronal networks: Effects of informa- 
tion transmission delay and rewiring probability, EPL 83 (2008) 
50008. 

[46] Q. Wang, M. Perc, Z. Duan, G. Chen, Synchronization transi- 
tions on scale-free neuronal networks due to finite information 
transmission delays, Phys. Rev. E 80 (2009) 026206. 

[47] A. Roxin, N. Brunei, D. Hansel, Role of delays in shaping spa- 
tiotemporal dynamics of neuronal activity in large networks, 
Phys. Rev. Lett. 94 (2005) 238103. 

[48] Y. Gong, Y. Hao, X. Lin, L. Wang, X. Ma, Influence of time 
delay and channel blocking on multiple coherence resonance in 
Hodgkin-Huxley neuron networks, BioSystems 106 (2011) 76. 

[49] Y. Hao, Y. Gong, X. Lin, Multiple resonances with time de- 
lays and enhancement by non-gaussian noise in newman-watts 
networks of Hodgkin-Huxley neurons, NeuroComputing 74 
(2011) 1748. 

[50] P. Balenzuela, J. Garcfa-Ojalvo, Role of chemical synapses in 
coupled neurons with noise, Phys. Rev. E 72 (2005) 021901. 

[51] N. Brunei, Dynamics of sparsely connected networks of ex- 
citatory and inhibitory spiking neurons, J. Comp. Neurosci. 8 
(2000) 183. 

[52] S. Lee, A. Neiman, S. Kim, Coherence resonance in a Hodgkin- 
Huxley neuron, Phys. Rev. E 57 (1998) 3292. 

[53] Q. Wang, X. Shi, G. Chen, Delay-induced synchronization tra- 
sition in small-world Hodgkin-Huxley neuronal networks with 
channel blocking, Discrete and Continous Dynamical Systems 
Series B 16(2011)607. 

[54] C. Koch, Biophysics of Computation: Information Processing 
in Single Neurons, Oxford University Press, Oxford, 1999. 

[55] Q. Wang, G. Chen, M. Perc, Synchronous bursts on scale-free 
neuronal networks with attractive and repulsive coupling, PLoS 
ONE6(2011)el5851. 

[56] V. Volman, I. Baruchi, E. Ben-Jacob, Manifestation of function- 
follow-form in cultured neuronal networks, Physical Biology 2 
(2005) 98. 

[57] C. Zhou, L. Zemanova, G. Zamora, C. C. Hilgetag, J. Kurths, 
Hierarchical organization unveiled by functional connectivity 
in complex brain networks, Phys. Rev. Lett. 97 (2006) 238103. 

[58] J. Gao, W. Tung, Y. Cao, J. Hu, Y. Qi, Power-law sensitivity to 
initial conditions in a time series with applications to epileptic 
seizure detection, Physica A 353 (2005) 613. 

[59] Y.-C. Lai, M. G. Frei, I. Osorio, L. Huang, Characterization 
of synchrony with applications to epileptic brain signals, Phys. 
Rev. Lett. 98 (2007) 108102. 

[60] V. Volman, M. Perc, M. Bazhenov, Gap junctions and epilep- 
tic seizures two sides of the same coin?, PLoS ONE 6 (2011) 
e20572. 



